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Abstract I outline a new model of particle acceleration in the current sheet sepa- 
rating the closed from the open field lines in the force-free model of pulsar mag- 
netospheres, based on reconnection at the light cylinder and "auroral" acceleration 
occurring in the return current channel that connects the light cylinder to the neutron 
star surface. I discuss recent studies of Pulsar Wind Nebulae, which find that pair 
outflow rates in excess of those predicted by existing theories of pair creation occur, 
and use those results to point out that dissipation of the magnetic field in a pulsar's 
wind upstream of the termination shock is restored to life as a viable model for the 
solution of the "a" problem as a consequence of the lower wind 4-velocity implied 
by the larger mass loading. [M 



1 Follow The Energy 

Rotation Powered Pulsars (RPPs) provide the first and most definitive example of 
compact astrophysical systems which draw the power for their observed emissions 
from the extraction of rotational energy from gravitationally bound objects through 
the action of macroscopic electromagnetic fields. They have motivated models for 
similar energy extraction from disks around other gravitating bodies, such as black 
holes (e.g. |Rees| ( |1984] >; |Begelman et al.\ ( |1984| l). As objects of study, the RPPs 



have a virtue lacking in the black hole systems: timing of the precisely measured 

Department of Astronomy, Department of Physics, Space Sciences Laboratory and Theoretical 

Astrophysics Center 

University of California, Berkeley 

e-mail: arons@berkeley.edu 

*With apologies to Charles Dickens 

* Collaborators, none of whom should be held responsible for the content of this paper: D. Alsop, 
E. Amato, D. Backer, P. Chang, N. Bucciantini, B. Gaensler, Y. Gallant, V. Kaspi, A.B. Langdon, 
C. Max, E. Quataert, A. Spitkovsky, M. Tavani, A. Timokhin 



1 



2 



Jonathan Arons 



pulse periods, uniquely interpretable as the rotation periods of the underlying neu- 
tron stars, provide measurements of the total energy budget free of all astrophysical 
uncertainties, other than the factor of ~ 2 uncertainty in neutron stars' moments of 
inertia, arising from the uncertainties in the equation of state of dense matter. The 
measured rate of rotational energy loss, 

E R = -IG*G* = 47tip, (1) 

tells us the total energy budget for these systems, without our having to understand 
anything about the photon emissions from these systems - which is both a blessing 
and a curse - a blessing, because in contrast to other relativstic astrophysical sys- 
tems, we know the energy budget, without having to unravel the partition between 
flow kinetic energy, large scale Poynting flux, thermal energy and radiative losses 
- a curse, since the energy loss is radiatively silent, thus supplying little informa- 
tion as to the details of the energy outflow, leaving the mechanics of the machine 
mysterious. 

Nevertheless, progress has been made. The Pulsar Wind Nebulae (PWNe) act as 



catch basins for the rotational energy lost. Observations of these systems Gaensler 
|& Slane| ( |2006| l, using radio (including millimeter), X-ray, gamma ray and occa- 
sionally infrared telescope^] have made clear that RPPs deliver their energy to the 
outside world in the form of highly relativistic, magnetized outflows - stellar winds 
that are exaggerated versions of the solar wind - which must be electromagnetically 
driven by the magnetic pressure of the wound up magnetic field. The strength of that 
field is estimated by using the theory of magnetic braking of the neutron stars' spin, 
which suggests 

G = -KG", G =2% IP, (2) 

P = rotation period, applied to the observed rotations periods and spindown rates 
P = -InG/G 1 . 

The earliest model applied vacuum electrodynamics to a rotating sphere endowed 
with a magnetic dipole moment fi centered at the stars' centers and tipped with re- 
spect to the rotation axis by an angle i. Tha t theory yiel d s the spindown luminosit y 



E R =KG 4 , K = (2 / '3) /I 2 sin 2 i/c? [e.g. |Pacini (1967i; Ostriker & Gunn <1969i] 



thus n — 3 in this model. Vacuum theory was motivated by the large gravitational 
forces at the surfaces, suggesting no plasma more than a meter or so above the 
star, but immediately Deutsch's much earlier observation [Deutsch]( |1955) ], made 
in the context of magnetic A stars, that in vacuum large electric fields parallel to 
B would overwhelm gravity and pull charged particles out from the star until the 
vacuum electric field would be altered, reducing E ■ B down to zero, was recov- 
ered (Goldreich & Julian 1969| l and extended with the suggestion that the charged 



particles would feed a curious charge separated wind, with a total electric current 



Most PWNe lie in the galactic plane, therefore are relatively inaccessible to optical techniques, 
and are even less accessible to UV telescopes. Near and far infrared observations are extremely 
useful in unraveling the physics of the relativistic outflows (Bucciantini et al. 2010 ), but have 
been much less in evidence than the high energy studies. 
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I = c<J> mag , <5 mag = total magnetospheric potential = ^Er/c. That wind could carry 
away the rotational energy of even the aligned rotator, in a Poynting flux dominated 
flow - the electromagnetic energy density would vastly exceed the kinetic energy 
density (and pressure) of the outflow, in the initially conceived model - the particle 
flux in that scheme is only c<f>/e = 2.3 x 10 30 (I45P15 / P 3 ) 1 / 2 elementary charges/s, 
7 45 = 7/10 45 cgs,Pi5 = the "Goldreich-Julian" current. 

The charge separated model has several really serious theoretical difficulties, but 
perhaps of greater importance is that the observations of young PWNe have made 




Period (s) 



Fig. 1 Observed RPP periods and period derivatives, from Kaspi • 2010 '. The line bounding the 
pulsar population corresponds to the magnetospheric voltage <P — 10 lz V, clearly marks a boundary 
beyond which pulsar emission is unlikely. Simple estimates of pair creation suggests this source 
of plasma should occur only for voltage larger than this value, which underpins the idea that pair 
creation is esential for radio emission |Sturrock 1 1971 )]. The change in slope of this "death line" 
al short period, small P indicates more sophistication in the pair creation physics, and/or in the 
association of pairs with radio emission, than is incorporated in the simplest models, a conclu- 
sion also apparent from the variety of quantitative problems with this widely accepted hypothesis 
dHibschman & Arons||2001||Medin & Lai||2010k. 



4 



Jonathan Arons 



clear that the particle outflows are many orders of magnitude larger than the 
Goldreich- Julian current, thus motivating MHD models, where E ■ B = is assumed 
from the start. Pair creation occurring somewhere within the magnetospheres (first 



suggested by (Sturrock 1971 1) is thought to be the origin of the required dense 
plasma, although a quantitative model that yields the high mass loss rates observed 
is still lacking [Bu cciantini et a/.| ( |2010| ) and references therein]. In recent years, 
solutions for the MHD structure of the magnetosphere in the appropriate force- 
free limit have been obtained (numerically - the attempts of analytically minded 
theorists to guess the answer to the rather formidable free boundary/eigenfunction 
problem posed by the steady state version of the model, uniformly failed over 30 
years of attempts), first for the aligned rotator [Conto poulos et q/.| ( |1999| ); |Gruzinov| 
d2005ll; IKomissarovl (|2006>; |Timokhin| (|20"0"6l>], then for the full 3D oblique rotator 



ISpitkovsky (2006); Kalapotharak os, C, & Contopoulos7Ll ( |2009[ )]. 

As far as the P, P diagram goes, the main result from these investigations is the 
innocuous looking conclusion in expression Q. 



K = k(l+sin 2 i) 



3, ,2 



k= IzbO.l. 



(3) 



Physically, the most important result has been the identification of the current sheets 
separating the closed from the open regions of the magnetosphere, extending into the 
wind beyond the magnetosphere, whose last closed flux surface ends just touching 
the light cylinder whose cylindrical radius is GJ = c/£2. That such a current sheet 
should be present has been suspected from the early days of RPP research [e.g., 
Michel] ( fl975| )]. The error in k reflects the uncertainties in the numerical treatment 
of the problem, many of which are associated with how the current sheet is repre- 
sented in the numerical schemes. Figure [2] shows a slice through the 3d force-free 
magnetosphere of the 60° rotator, from Spitkovsky's (2006) results. Observationally, 
expression Q shows that the vacuum rotator's "braking index" n — 3 is preserved in 
full force-free MHD, which contradicts the observed values in the small number of 
stars where n has been determined (Livingstone et al. , 2007), a contradiction which 
has led to a variety of suggestions ranging from evolution of the magnetic moment 
fl or the obliquity i (Blandford & Romani 1988) to effects of reconnection on the 
rate of conversion of open magnetic flux to closed ( |Contopoulos & S pitkovsky, 



2006), as the star spins down and the closed zone expands at the expense of the 



amount of open magnetic flux. That reconnection might affect the braking index 
is readily derived from the fact that the torque really depends on the magnitude of 
the open magnetic flux. The amount of open flux depends on the size of the closed 
zone, which ends at Ry. If Ry/Rl < 1, the torque increases because of more open 
field lines and larger Poynting flux than is the case for a magnetosphere closing at 
r = Ri. Bucciantini et al. ( 2006[ ) show that the braking index is 
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Fig. 2 Field lines and current density of the oblique force-free rotator, i = 60°, from Spitkovsky 
(2006), seen in a cut in the Q , fl plane. The last closed field lines end at a Y-point (a Y-line, in 3D) 
at distance Ry from the neutron star. The current sheet encloses and separates the closed from the 
open field regions of the magnetosphere, and the separate branches merge in the wind zone, where 
the folded sheet continues to separate the oppositely directed fields of the striped wind. The arrows 
and "Gap" labels locate sites where vacuum gaps have been postulated, in test particle models of 
accelerators that lead to gamma ray emission and pair creation. None of these gaps appear in a 
current sheet accelerator based on the force-free magnetosphere model. 



base of the current sheet in the relativistic wind of the aligned rotator. The specula- 
tion is that spindown slightly biases these reconnection events so that on the much 
longer spindown timescale, the net open flux slowly converts to closed. These cur- 
rent fluctuations might be associated with the timing "noise" ( Arons 1981a JCheng, 



et al. 



1987 ) identified long ago with torque fluctuations (e.g., Helfand et aOi T980F Scott 



( 2003) >), although recent analysis of longer data sets fLyne et al. |2010 1 in long 



period pulsars has called the noise interpretation into question. 

In itself, the force-free model does not provide mechanisms for photon emission. 
But it has a variety of implications, which are slowly being addressed. 



• The model specifies the polar flux tube size and shape - it is noncircular with a po- 
lar cap center displaced from the magnetic axis, even when the magnetic field is 
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the simplest, that of a star centered, point dipole (Bai & Spitkovsky 20101. This 
has consequences for radio polarization structure, and for polar cap areas and 
dipole offsets inferred from soft X-ray emission from polar caps (Bogda nov et\ 
al. 2007|l, thought to be h eated by magnetospheric particle bombardment ( Arons , 
1981b[|Zavlin & Pavlov[ [1998] parding & Muslimovl [2002] ). These theoretical 
improvements of the polar cap model have yet to be noticed and incorporated in 
phenomenological models of the observations used by data analysts; such incor- 
poration might yield interesting tests of the force-free model. 



Dip: magnetic field 
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Fig. 3 Plasmoids formed at the base of the current sheet of the relativistic aligned rotator, where 
the sheet crosses the light cylinder, from a relativistic MHD simulation of a newly born neutron 
star's magnetosphere (Bucciantini et al. , 2006). The dissipation that allows reconnection to occur 
is numerical. They move out radially at the local Alfven speed Va ~ c, and recur on the magneto- 
spheric Alfven transit time ~ P jiz 



• The force free model quantitatively specifies the return currents required to pre- 
vent the star from charging up, as the polar flow extracts charge from the star. 
The results for the oblique rotator are partly in accord with long held expecta- 
tions, that return current exists in a thin ("auroral") sheet bounding the polar flux 
tube {e.g. (Goldreich & Julian 1969 |Michel| [l975| l), consistent with the open 

1969| l - current closure occurs far away, in 



circuited model ( Goldreich & . 



ulian 



the nebula/interstellar medium beyond the wind termination shock or perhaps in 
the outer wind, plus qualitatively new features: a) part of the return current sur- 
rounding the polar flux tube is spatially distributed, even in the aligned rotator; 
b) in the oblique rotator, this part of the current system is not a return current at 
all, but couple the two polar regions together ( |Bai & S pitkovsky] |2010| l - in the 
orthogonal rotator (z = 90°), the auroral component of the current is entirely in 
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the polar coupling flow, with the volume current out of each half of the polar cap 
having equal amount and opposite sign, also consistent with early expectations 



(e.g. Scharlemann et. al. ( 1978 1) - the orthogonal rotator automatically balances 
it's charge loss. The radiative consequences of these features are as yet unex- 
plored - for example, the spatially distributed part of he return current might be 
a good candidate for the site of "conal" component of pulsar radio emission, an 
idea which requires non-force modeling of the current flow and identification of a 
workable emission process within that current flow modej^Jbefore one can relate 
the theoretical force free current distributions to the observations in a testable 
manner (although easier kinematic comparisons are certainly possible). 
The location of the return current layer having been determined, the hypothesis 
that the return current layer is the site of the beamed particle accelerator that 
gives rise to the pulsed gamma rays observed by the FERMI and earlier orbiting 



gamma ray telescopes (see|Ray & SazParkinson (20101 in these proceedings for 



a recent review) can now be investigated in the context of a self consistent mag- 
netopspheric structure that allows a quantitative evaluation of the beaming char- 
acteristics implied by the radiating current sheet concept - see ( |Bai & Spitkovskyl 
201 0| l for a kinematical study of the radiating current sheet idea. 



Making progress on a physical model for radiation from the current layers can 
most expeditiously take advantage of the facts that a) pulsed gamma ray emission, 
when observed, is the largest photon output from rotation powered pulsars, but b) 
generally has less luminosity than the spin-down luminosity of these stars. Figure 



2 Some earlier ideas on this subject relating to field aligned acceleration and gamma ray emission 
can be found in |Arons|[l983b||Gruzinov||2008| , for example. 
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[4] illustrates this fact, which summarizes the results from the LAT instrument as of 
Spring 2010: 



35 



l 34 

i 



33 



4V 



t 



Laccel/Lspindown << I 



34 



35 36 37 

Log[E (ergs" 1 )] 



38 



39 



Fig. 4 Ratio of observed gamma ray luminosity L gamma to observed spindown power E = Z. sp i n down 
for LAT pulsars {Abdo et q/.||2010) . If the acceleration of the gamma ray emitting particles is ra- 
diation reaction limited (generally true in models that assume curvature radiation as the emission 
mechanism, as is the case in most of the "gap" models). The trend at higher luminosities roughly 

follows the efficiency Lg amma /L spm down x (Vc expected if some mechanism limits the acceler- 

ftjc, [e.g 



ating voltage A<t> to a fixed fraction of the magnetospheric voltage <P mi g 
( 1996); Harding & Muslimov {2002 )], a limitation plausibly attributed to pair creation. The large 
dispersion in Figure |4j comes primarily from the uncertain distances, and also to the uncertain 
beaming correction required to go from the fraction of the sky illuminated by the beam to the total 
emission that would be detected by an observer covering An ster - the LAT team assumed uniform 
phase-averaged beaming across the sky (1 sterradian). 



Arons 



Thus the energy invested in particle acceleration is a small fraction of the enegy 
stored in Poynting fluxes, for radiation reaction limited acceleration, so that the force 
free model can be considered as a good zeroth order magnetospheric description. 

The voltage limitation prediction L gamma oc <J> mag oc y / £ is the same as L gamma °c 
"particle current" = particle flux = Goldreich-Julian flux in unidirectional beam 
models of the polar electric current flow (Harding 1981| l only if the accelerator 
carries a fixed fraction of the total electric current that enters into the spin-down 
torque, independent of P,P. In the traditional slot or outer gap models, pair creation 
establishes the limitation of the accelerator to a thin sheet either in the outer mag- 
netosphere ( |Cheng, Ho & Ruderman| [T986 ) or back in the polar cap ( |Arons| [T983a) 
- in such models, the gap width w* of a model which successfully reproduces the 
sharply peaked light curves, projected onto the neutron star following the field the 
poloidal field lines, is necessarily small compared to the polar cap size, and varies 
with the pulsars' magnetic moment and spin parameters. This is not an issue for a 
model based on the currents flowing in the return current layer, which necessarily 
carries the whole magnetospheric return current and include a substantial fraction 
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in a thin current sheet, for most obliquities of the magnetic moment with respect to 
the rotation axis. 

That magnetospheric current sheets, with particle densities within the sheets high 
compared to the Goldreich- Julian value, can sustain large parallel electric fields is 
well known in planetary magnetopsheres - such sheets are the accelerator sites of the 
particle beams that stimulate the aurora observed in the upper atmospheres of the 
Earth, Jupiter and Saturn, for exampl^An elementary illustration of this possibility 
comes from considering the inertial term in the generalized Ohm's for the electric 
field parallel to B, which can be written as, in the relativistic case, as 

E H = -2 + l • V- [7(JV + VJ)]} oc (5) 

(Op \ dt B "} 11 ZWientPfi 

with my the particles' relativistic mass, n their density, / the total current set by the 
force free magnetosphere, 4 curren t is the thickness of the current carrying channel, 
and Pb is the radius of curvature of the magnetic field. Thus inertia of the current 
carriers can act as a effective resistance in these high inductance systems, which 
establish the currents electromagnetically, forcing a parallel electric field to appear 
in the presence of any parallel (to B) load - pressure and radiation reaction itself are 
other effects which can serve as loads, with pressure especially important in the dif- 
fusion region around the singular line shown in figure[5] the fluid velocity. Equation 
|5]) is most useful in the corotating frame, and when the current is due to relative mo- 
tion between the species (electrons and positrons, and in some circumstances heavy 
ions) which is slow compared to the bulk fluid velocity v. In the current circum- 
stance, it turns out that the current is better described as counter-streaming beams 
- in that case, describing the beams as separate fluids is more appropriate, and they 
can form the total plasma density in the current flow channel, rather than being a 
low density component in a much denser plasma. 

Expression <|5j does make clear that acceleration is prone to maximize when the 
the relativistic mass is high and the current density is high (large I/A current ). A cun:ent 
is generally microscopic, expected to be on the order of c/(O p , and is established 
by the dynamics of the the singular region where the closed zone ends, illustrated 
(in the cartoon approximation) in Figure [5] The capture rate of pair plasma into the 
diffusion region is 

Vrin ^ 2l D A L fe c Cfinag ,„ 

- ,v diffusion± tfT fl ± Z, ' W 

L Pwind " 

with K± the multiplicity (multiplier of a fiducial Goldreich- Julian outflow rate 
c<5 ma g/e) that gives the number of pairs in the total outflow, p\ ec the reconnection 
speed in units of c and j3 w ; n( i the polar wind outflow velocity (set equal to unity). Al 
is the thickness of the current channel at the light cylinder, assumed equal to the half 
height of the diffusion region, and Id is the length of the diffusion region. Pressure 



For a review of such phenomena, see Paschmann et al. • 2002 i. Of particular significance to 



pulsars is the fact that the field aligned currents that power narrow auroral arcs consist of pre- 
cipitating electron beams launched from the reconnection region in the distant magnetotail and 
counterstreaming ions launched from the planetary atmosphere. 
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in the diffusion region expels the captured pairs from the diffusion region, outwards 
along the current sheet in the wind and inwards along the auroral current channels. 
Expression (|6|l assumes the plasma flux in the wind has a gradient across B, as is 
likely since the accelerating electric field in the polar cap that leads to the pair cre- 
ation that feeds the wind is small near the cap edge. The pressure supported electric 
field provides the accelerator which sorts the particles in the diffusion region into a 
precipitating beam (electrons in the geometry shown in Figure [5]l and an oppositely 
charge beam traveling outwards in the wind current sheet - the charge signs of the 
beams are as required by the global electrodynamics. The channel thickness is al- 
most certainly comparable to the skin depth in the pairs. Taking the plasma gradient 
into account leads to the lower limit to the channel thickness 



co p ±(5 = c/co p ±(5) \2K ± e<P miig 

where 8 is the distance across B from the formal current sheet location, and a pre- 
cipitating beam flux in the return current channel 

Fv4^W; (8) 
R L 2nR{ \r J 

in the gradient case, Al drops out of the precipitating flux evaluation. y± is a mea- 
sure of the four velocity and of he comparable four- velocity dispersion of the po- 
lar plasma flow emerging from the inner magnetosphere, predicted by pair cascade 
models to be on the order of 10 2 . Numerically, ^ yields a very small value, on the 
order of meters to hundreds of meters, the specific value depending on 4>mag and 

The kinetics of relativistic reconnection being a largely untrodden subject, Id is 
more or less unknown - it could be as small as Ai itself (Petscheck style reconnec- 
tion), or as large as appears in the numerical dissipation driven reconnection ob- 
served in the force free simulations, Id ~ O.IRl (A. Spitkovsky, personal communi- 
cation). In the terrestrial magnetosphere, satellite observations suggest the diffusion 
region length is intermediate between the ion skin depth and the macroscopic scales. 
For the discussion here (and in the more detailed report in preparation), treating Id 
as a parameter to be constrained by model comparisons to observations appears to 
be the wisest strategy. In analogy to observed non-elativistic reconnection, one ex- 
pects j3 lec ~ 0.1vAifven/ c = 0.1, a value supported by the few PIC simulations of 
relativistic reconnection [e.g., Zenitani & Hoshino (20(37]!]. 

The precipitating particles form the "hanging charge clouds" invoked by |Goldre^| 
|ich & Julian| (fl969) to cause the electrostatic extraction of return current from the 
star's atmosphere, which happens if « v (^*) = Fy(R*)/c ^> ncj(R*), or , from d8j, 
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Fig. 5 Upper Panel: Electric current structure of the oblique force-free magnetosphere, inclination 
= 60°, from Spitkovsky ( 2006 1. The current sheet, indicated by the darker color, bounds the closed 
zone. The closed zone ends at a singular Y-line (in the ideal force-free electrodynamics approxi- 
mation) at the light cylinder distance from the neutron star. Lower Panel, left: Current structure at 
the polar cap, illustrated for acute angle between the magnetic moment and the angular velocity 
(electron polar current, Q ■ fl > 0), The return current in the current sheet consists of a precipitat- 
ing electron beam, launched from the diffusion region around the Y-line and possibly augmented 
by high altitude pair creation within the current sheet, with the charges in the precipitating beam 
extracted by reconnection flow from the pair plasma flowing from the polar cap into the wind, 
plus a counterstreaming ion beam extracted electrostatically from the stellar atmosphere. In obtuse 
geometry (ion polar current, Q. ■ /! < 0), positrons precipitate from the Y-line and counterstreaming 
electrons are extracted from the atmosphere. For clarity, the part of the return current not contained 
in the current sheet is omitted, even though this part of the current system is of increasing signifi- 
cance as i — > 90°. Lower Panel, right: Possible structure of the Y-line region, with the termination 
of the closed zone to the left and the merger of the winds from the opposite polar caps to the right. 
The "guide field" also reverses across the mid-plane of the flow, along with the poloidal open 
field. As reconnection occurs, some of the ouflowing plasma (speed cj3, /3 « 1) deflects toward the 
singular, unmagnetized "diffusion" region around the Y-line with speed v rec ~ O.h'Aifven = 0.1c. 
The figure represents a steady (in the co-rotating frame) flow model - in reality, the reconnection 
is likely to be bursty, as in Figure^] with formation of sporadic X-lines. 
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if 

^f^K ± »L (9) 

K L Pwind 

The condition in (|9]l is satisfied if lo is large compared to the skin depth, but still 
small compared to the numerical dissipation determined length observed in the force 
free simulations. It can be shown that the consequent minimum total potential drop 
in the twin beam channel is 

1 /3 

A ~ 1 ^ R * f m ± c 2 r± /3 w ind\ Id 

8 *R L \2e<P mag K ± ^ / 4 L 



1^ R* /jSwindV 73 jg 



-O^knagTH -75 ) ~ cos '- (10) 



If Id/Rl really is as large as O.l, rather than comparable to the minimum scale Al, 
the accelerating potential in the return current channel is more than large enough 
(around 10 13 Volts, in the Crab and Vela pulsars) to drive curvature gamma ray 
emission from the beams, which limits the particle energies by radiation reaction. 
Such energies are high enough to lead to pair creation, which may limit the acceler- 
ation, although as is clear from expression Q, parallel potential drops in the current 
carrying region can be sustained even if the plasma is dense. 

If one models the gamma ray emission as being directly from the beams, one 
needs to take account of the electron beam going inwards, in acute geometry, where 
the upward beam is ions (generically, protons) - for obtuse geometry, the outbound 
beam is electrons, with an inbound positron beam. Pair creation (through y— y in- 
teractions with softer photons from the star), when it is important, makes the obtuse 
and acute geometries equivalent. 

The radiation might be synchrotron emission, especially the nonthermal X-ray 
output. The counterstreaming beams are electromagnetically two stream and shear 
unstable. Since (Qpbeam can be comparable to the relativistic cyclotron frequency 
in the outer magnetosphere, the growing waves can excite finite Larmor gyration 
of the particles in the current channel, thus producing incoherent emission through 
hard X-rays - gamma rays are possible, under some circumstances. These X-rays are 
an alternative to soft photons from the star, as targets for J— J pair production. If the 
lower frequency waves can escape the plasma, they are a direct source of coherent 
emission, perhaps of interest to modeling giant radio pulses, which appear to come 
from the outer magnetosphere. 

A more detailed description of the model will appear elsewhere. 



2 Follow the Mass 



Pulsar Wind Nebulae (PWNe) demonstrate that pulsars loss rest mass at a rate large 
compared to the fiducial electrodynamic particle loss rate c<J> mag /e [for a review, 
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see |Gaensler & Slane| p006 )]. The only known explanation is pair creation in the 
pulsars' magnetospheres. In many nebulae, the X-ray emitting particles rapidly lose 
energy to synchrotron radiation. Then the nebulae are particle and energy calorime- 
ters, allowing direct inference of the pair multiplicity in the wind, of TeV to PeV 
pairs. The measured injection rates, up to ~ 10 385 pairs/s, compare well to the pre- 
dictions of existing pair creation models [e.g. Hibschman & Arons (2001 )], yielding 
multiplicities up to ~ 10 4 . However, PWNe are also radio synchrotron emitters, ra- 
diation that samples much lower energy populations (100 MeV to 10 GeV), whose 
radiative efficiency is much less than their X-ray emitting cousins. The result is a 
much larger population of pairs, whose radiative lifetime exceeds that of the nebu- 
lae. The most efficient hypothesis is that these particles come from the embedded 
pulsars also, an idea supported by spectral continuity and by exotica, such as the 
observation of radio "wisps" near the Crab pulsars (Bietenholz et al. 2001 1. Apply- 
ing simple evolutionary models allows one to infer time averaged injection rates. 
Recent evaluations by Bucciantini et al. ( 20 1 0) > and ( |Slane et ~al. 2010 1 yield lower 
limits for multiplicities K± all in excess of 10 s and upper limits for wind 4-velocities 
Iwind = Er/Mc 2 = e<t> mag /2K±m±c 2 all less than 10 5 in a number of nebulae. The 
data are the best for the younger systems, although even for these, the lack of far 
infrared data inhibits the analysis. 
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The inferred multiplicity excesses are a puzzle for pair creation theory [e.g. |Hi- 
|bschman & Arons] ( |2001| l; |Timokhin| ( |2010| l], perhaps resolvable by appealing to 
magnetic anomalies near the neutron stars' surfaces, the simplest being an offset 
of the dipole center from the stellar center, which strengthens the magnetic field 



at one pole (Arons 19981. The increase this gives to the magnetic opacity can be 
greatly enhanced if the magnetic axis is also tipped with respect to the radial di- 
rection, since then gravitational bending of photon orbits with respect to the B field 
direction much increases the magnetic opacity for pair creation. Such phenomeno- 
logical modifications of the low altitude magnetic field must respect the observation 
that radio beaming morphology is consistent with the magnetic field being that of a 



star centered dipole quite close to the star ( |Rankin 1990; Kram eref q/.| [l998 ). This 
problem warrants quantitative investigation. 

The large inferred multiplicities imply the wind 4- velocities r^j n d to be small 
compared to the much quoted value of 10 6 inferred by Kennel & Coroniti ( 1984 1 



in their model of the Crab Nebula's optical and harder emission. The large mass 
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loading and the inferred low wind four velocity has a large impact on the much sto- 
ried "a problem" of pulsar winds. In ideal MHD, the ratio a of magnetic energy to 
kinetic energy in the wind is conserved outside the fast magnetosonic radius (since 
for a cold flow the wind does not substantially accelerate outside this surface) and 
is large - even with the increased mass loading found from recent nebular studies, 
a is always well in excess of several hundred. Nevertheless the wind behaves at its 
termination shock as if a is small - MHD models of the nebulae suggest a at th 



termination shock is on the order of 0.02 in the Crab Nebula [e.g DelZanna et al. 
( 2004|>] and similar values are plausible in other systems 



Coroniti ( 1990) suggested that because the wind of an oblique rotator has the 



magnetically striped structure shown in Figure [6] magnetic dissipation of the corru 
gated B field, generically of a resistive nature propelled by instabilities of the current 
flow in the current sheet separating the stripe^] might destroy the magnetic field of 
the wind interior to the termination shock, thus converting a high a flow into an ef- 
fectively unmagnetized plasma. If the current sheets separating the magnetic stripes 
are to merge with a speed v s < c as measured in the proper frame of the flow, be- 
fore they reach the termination shock located at distance Rjs from the neutron star, 
the merger time in the PWN frame r merge = ^r^: ind (Ri/v s ) must be less than the 
flow time to the termination shock Rjg/c, therefore -T^j n( j 

5 x 10 4 y/v s /c (Crab) must be satisfied if Coroniti's model is to be viable. This in- 
equality is satisfied for multiplicities above 10 5 , which does appear to be the case 
for the young PWNe recently analyzed. In this low i^vi nt j, high pair multiplicity en- 
vironment, complete destruction of the striped component of the wind's magnetic 
field is possible ( Arons 2008 ), contrary to the conclusion reached by (|Lyubarsky & 



Kirk 2001 1, but in accord with the conclusions of ( |Kirk & Skjaerassen, 2003 i, al- 



though for somewhat different reasons (drift kink instability dominate over tearing 
in pair plasmas). 

Possible mechanisms that can lead to the necessary dissipation are drift-kink 



instability of the current sheet, considered as if it were a flat sheet (Zenitani & 



|Hoshino||2007| l and an interesting Weibel-like instability due to interaction between 
the sheets (Arons , 2008 ), an effect strongest in the equatorial sector where the folded 
sheet appears locally as neighboring flat sheets with antiparallel current flow in the 
latitude direction. Most interesting from the astronomical perspective, the dissipa- 
tion of the current sheets in the inner wind (r <C Rts> while responsible for only a 
small fraction of the magnetic destruction, may have sufficient luminosity to allow 



This wrinkled current sheet, frozen into the wind, is the second current sheet of this paper's title, 
but really it is the continuation of the sheet separating the closed and open zones interior to the 
light cylinder, as is apparent in Figure[5] 




Fig. 6 a) Magnetic Geometry of a Force-Free Rotator for r < 2Ri, for j = 60°, from Spitkovsky 
(2006). The rapid transition to inclined split monopole field geometry for r > Rl is apparent, b) 
Geometry of the current sheet from the split monopole model for i = 60°, r > Rl. For clarity, only 
one of the two spirally wound current sheets is shown. As i — > 90°, the sheets almost completely 
enclose the star; for r^>Ri, the spirals are tightly wrapped (B r <C B^) and the current sheet surfaces 
closely approximate nested spheres, c) One sheet for = 30°, shown for clarity, d) Meridional cross 
section of the current sheet for / = 60° . e) Equatorial cross section snapshot of the current sheet, 
showing the two arm spiral form. The arrows show the local directions of the magnetic field; 
the dots and crosses show the direction of the current flow. Panels b)-e) were constructed using 
Bogovalov's ( Bogovalov] ! 1999} analytic model of the asymptotic wind, f) Current sheet from a 2D 
PIC simulation of the inner wind, from unpublished work by Spitkovsky (used by permission), 
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